Velocity distributions in dissipative granular gases 



J. S. van Zon 1 ' 2 and F. C. MacKintosh 1 - 2 

1 Division of Physics and Astronomy, Vrije Universiteit, 1081 HV Amsterdam, The Netherlands 
2 Kavli Institute for Theoretical Physics, University of California, Santa Barbara, CA 93106 

(Dated: February 2, 2008) 

Motivated by recent experiments reporting non-Gaussian velocity distributions in driven dilute 
granular materials, we study by numerical simulation the properties of 2D inelastic gases. We 
find theoretically that the form of the observed velocity distribution is governed primarily by the 
coefficient of restitution n and q = Nh/Nc, the ratio between the average number of heatings 
and the average number of collisions in the gas. The differences in distributions we find between 
uniform and boundary heating can then be understood as different limits of q, for q 1 and q < 1 
respectively. 
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Granular materials consisting of macroscopic particles 
or grains can exhibit behavior reminiscent of conventional 
phases of matter. Sand, for instance, can flow like a liquid 
under some conditions. Dilute granular systems, or gases, 
have been extensively studied both experimentally and 
theoretically, in large part as simple model systems ex- 
hibiting nonequilibrium and dissipative behavior. These 
systems are intrinsically dissipative and out of equilib- 
rium, even though it is tempting to apply such equilib- 
rium notions as temperature. Since the collisions in such 
a gas are inelastic, a gas-like steady state is achieved only 
with a constant drive, or input of energy. Otherwise, all 
motion ceases after only a finite time 0,0- In principle, 
it is possible to drive the system uniformly throughout 
the container, with every particle in contact with a heat 
source all the time (uniform heating). This has been done 
in simulations 0, Q and is assumed in analytic theories 
0. In experiments, however, one usually drives a granu- 
lar gas by shaking or vibrating the walls of the container. 
Such boundary heating means that the e nerg y is inserted 
in a spatially inhomogeneous way 0, Ej- As a con- 
sequence, the gas will develop a gradient in density and 
mean kinetic energy @ . Even for uniform heating, how- 
ever, significant deviations from equilibrium gases, e.g., 
in density correlations, are observed 

One of the most fundamental aspects of molecular 
gas es is the Maxwell-Boltzmann velocity distribution 
A very general and striking feature of driven dissi- 
pative gases, however, is the apparently strong deviation 
from this classical behavior. Such velocity distributions 
have been the subject of numerous recent experimental 
studies 0, H, U El- an( i characteristic velocity distri- 
butions of the form P(v) — Cexp[— /3(v/a) a ] have been 



observed, where a 



')2 is often called the granular 



temperature, in analogy with equilibrium gases. Experi- 
ments of Rouyer and Menon |8|, in particular, have sug- 
gested that a universal exponent of a = 1.5 occurs over 
a wide range of experimental parameters. This observa- 
tion was particularly intriguing, as Van Noije and Ernst 
0] predicted a high- velocity tail with an exponent 1.5 us- 
ing kinetic theory. Here, we demonstrate instead that a 
family of distributions showing apparent exponents cov- 



ering a wide range of values a < 2 is expected, depending 
on both material and experimental conditions. Further- 
more, we show that the velocity distribution is governed 
primarily by the relative importance of collisions to heat- 
ing, i.e., the way in which energy flows through the sys- 
tem of particles. Specifically, we introduce a new pa- 
rameter q — Nh/Nc, which measures the ratio between 
numbers of heating events and collisions experienced by 
a typical particle. These theoretical observations can ex- 
plain both the observed non-Gaussian behavior, as well 
as the ambiguities in the experimental and theoretical 
literature on dissipative gases to date. We also demon- 
strate that the behavior of the velocity distributions can 
be captured quantitatively by a simple model that takes 
only 77 and q into account, with no spatial degrees of 
freedom. 

Even though the distribution of velocities is such a ba- 
sic property of gases, results for velocity distributions in 
dissipative systems have remained ambiguous. Both in 
simulation and experiment, different setups and driving 
mechanisms usually give different behavior of the velocity 
distribution. For a setup where particles on a horizontal 
plate were driven in the vertical direction, Olafsen and 
Urbach found a crossover from exponential to Gaus- 
sian distributions as the amplitude of the driving was 
increased. The result of Rouyer and Menon Q was ob- 
tained for a different configuration where particles were 
confined between two vertical plates and driven in the 
vertical direction. Blair and Kudrolli 01 use a differ- 
ent setup where particles move along an inclined plane. 
Friction with the plane during collisions reduces the ef- 
fective coefficient of restitution to n « 0.5, in contrast 
to r] = 0.93 in the setup of Rouyer and Menon. They 
find the distribution with exponent a — 1.5 only in the 
very dilute case. Otherwise, the distributions deviate 
strongly from both Gaussian and the distribution ob- 
tained by Rouyer and Menon. For the denser case Blair 
and Kudrolli find a distribution that has one exponent 
for low velocities and crosses over to another for high ve- 
locities. This crossover is also observed in simulation by 
Moon, et al. ||. 

Thus, it remains unclear what the analogue (if any) 
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of the Maxwell-Boltzmann distribution is for dissipative 
gases, let alone what the origin of the non-Gaussian be- 
havior is. Puglisi, et al. have suggested that the de- 
viations are caused by the spatial correlations in the gas. 
They propose that for a region of uniform density the ve- 
locity distribution of the gas actually is Gaussian, with a 
density dependent width. The spatial correlations cause 
density fluctuations and they claim the non-Gaussian dis- 
tributions arise as an average over the velocity distribu- 
tions over these regions of different density. Experiments 
0, ^2 show that this is not the case. Even within win- 
dows of uniform density the velocity distributions remain 
non-Gaussian. Our results below confirm this. 

We study velocity distributions using an event-driven 
simulation of N particles of radius r moving in a two- 
dimensional box. Particles gain energy by heating and 
lose energy through inelastic collisions. When two parti- 
cles i and j collide their final velocities depend on their 
initial velocities in the following way: 

Vj = Vj — (Vi • Ty - Vj • ry ) Tij , ( 1 ) 

where rj is the coefficient of restitution and fy is the unit 
vector connecting the centers of particles i and j. 

When heating uniformly, we add every time step At 
a random contribution to the velocity of all the particles 
and we use periodic boundary conditions to simulate bulk 
behavior. Note that this is significantly different from the 
spatially homogeneous heating used in experiments [l^ . 
as here in uniform heating the driving is uncorrelated in 
space and time. When heating through the boundary, 
particles are confined in a box with radius R = 1. The 
particles receive a kick upon collision with the bound- 
ary. We assume that the collision between particles and 
the boundary is elastic and we add a random contribu- 
tion to the velocity in the direction perpendicular to the 
boundary. For more details, we refer to [l5|. We allow 
the system to reach steady state before taking data. For 
both uniform and boundary heating, data is taken peri- 
odically every At. 

One of the first striking differences between uniform 
heating and boundary heating is clustering. When heat- 
ing through the boundary, a stable liquid-like cluster sur- 
rounded by a hot gaseous state will form for low coeffi- 
cients of restitution r\ or high area fraction <f>. This occurs 
as particles are compressed in the center of the box by 
particles moving in from the boundary. The increase of 
density leads to collapse and a stable cluster is formed. A 
typical example is shown in Fig.Q] These clusters do not 
occur in our simulations with uniform heating. This is 
because particles are heated all the time, which prevents 
the collapse to a cluster. 

Velocity distributions obtained for uniform and bound- 
ary heating are shown in Fig. [21 When heating through 
the boundary, the gas develops a gradient in area frac- 
tion and granular temperature. In this case, we record 
the velocity distribution in a ring of radius 0.4 < r < 0.6 
around the center of the box, where the granular tem- 



perature is approximately constant. The velocity distri- 
butions in the different rings varies only slightly. For 
uniform heating, we collect data everywhere in the box. 




FIG. 1: Snapshot of a clustered state state for N = 350, 
4> = 0.05 and 77 = 0.6. Particles are indicated by circle, while 
the lines show the direction and magnitude of the velocity. 
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FIG. 2: Velocity distributions for N = 350 and cfi = 0.02. 
Shown are both results for uniform heating with 77 = 0.8(o), 
77 = 0.1(D) and results for boundary heating with r\ = 
0.9(0), V = 0.6(V) and 77 = 0.4(*). (a) P(v x /a x ) (b) 
— ln{—\n[P(v x /<Jx)]} versus ln(v x /<Jx)- A Gaussian is shown 
as a solid line. 

Fig. HI shows clear qualitative differences between uni- 
form and boundary heating. For uniform heating the ve- 
locity distribution is close to Gaussian for a large range 
of 77. In contrast, for boundary heating the distribution 
is only Gaussian in the nearly elastic case r\ = 0.9. As 
the coefficient of restitution is lowered, a crossover devel- 
ops where the exponent changes from a.\ ks 2 to a lower 
value a2- The exponent 012 becomes smaller as dissipa- 
tion is increased (for smaller 77 or higher <p) and we find 
any value in the range 0.7 < 0:2 < 2. As 02 decreases 
it becomes increasingly difficult to describe the high ve- 
locity tail with a single exponent ai. This crossover in 
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exponent is observed in boundary heating for all values 
of <fr and N. 

For a certain range of parameters we also find high- 
velocity tails with exponent 012 = 1.5. For their experi- 
ments Rouyer and Menon used N particles with ri_~ 0.9, 
where 100 < N < 500 and 0.05 < <f> < 0.25 Q. In 
Fig. we plotted the velocity distribution for rj = 0.9, 
<j> = 0.05 and several values of N. We also show the fit 
with a = 1.52 as made in This line clearly coin- 
cides with the high-velocity tail of the velocity distribu- 
tion found by the present simulation. This suggests that 
instead of a universal distribution with a = 1.5, they 
might have observed a part of a more complex velocity 
distribution, with more than one apparent exponent. 

The main difference between uniform and boundary 
heating is that in the first case heating takes place ho- 
mogeneously throughout the box, whereas in the latter 
case energy is injected inhomogeneously at the bound- 
aries. This is not the direct cause for the difference 
in velocity distributions. When heating homogeneously, 
one can go from a Gaussian distribution to one with a 
crossover by increasing the time between heatings |l5j 
allowing the average number of heatings per unit time 
to become smaller than the average number of collisions. 
The reverse is also true. When heating inhomogeneously 
through a boundary, one finds Gaussian distributions in 
very dilute systems, when the particles on average collide 
more often with the boundary than with other particles. 
This suggests that in our system the shape of the velocity 
distribution is not a function of all parameters 77, cf>, N 
and the details of the energy injection, but only of 77 and 
q = Nh/Nc, the ratio between the average number of 
heatings Njj and the average number of collisions Nc- 

This idea can easily be tested in boundary heating. 
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FIG. 3: (a) - ln{- ln[P(v x /a x )]} versus \n(v x /o x ) for N = 
350, <t> = 0.05 and r\ = 0.9 (o), N = 500, 4> = 0.05 and rj = 0.9 
(□), N = 350, <f> = 0.05 and 77 = 0.8 (o), N = 350, <j> = 0.25 
and rj = 0.9 (A). The solid lines correspond to the fit as made 
by Rouyer and Menon and has an exponent a — 1.52. The 
range of the solid lines corresponds to half the range used by 
Rouyer and Menon in their fit, but contains about 80% of 
their data points. 



When increasing the number of particles N or the area 
fraction <p, the average number of collisions increases. 
One can show in a mean field approximation that q ~ 
(Ncj))- 1 / 2 . The average distance a particle travels be- 
tween collisions is given by l co u ~ 1/0. For a box 
of area A the average distance between boundaries is 
given by ly ieat ~ A 1 / 2 ~ (7V/0) 1 / 2 . Finally, we know 
Nh/Nc ~ hoii/lheat- Our simulation obeys this approx- 
imation very well. In Fig. 0] we show velocity distribu- 
tions for 77 = 0.8 and different combinations of N and c/>. 
We measure the heating-dissipation ratio q in the sim- 
ulation and show velocity distributions with the same q 
on top of each other. For q = 1.3 and q = 0.13 we find 
excellent collapse for different iV and 0, even when we 
scale the system by a factor 8. For q = 0.013, where 
spatial correlations become very strong, we still find rea- 
sonable collapse. As we increase q we observe the usual 
pattern, where a crossover appears in a distribution that 
was initially close to a Gaussian. 
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FIG. 4: Velocity distributions for different values of the 
heating-dissipation rate q. Distributions with the same q 
are shown on top of each other. (A) q = 1.3 and we show 
N = 100 and <j> — 1 ■ 10 _3 (o), N = 200 and <j> = 5 • 10~ 4 (D), 
N = 800 and <j> = 1.25 ■ 10~ 4 (0). (B) q = 0.13 and we 
show N = 100 and <j> = 0.08(o), N = 200 and cj> = 0.04(D), 
N = 400 and (j> = 0.02(0). (C) q = 0.013 and we show 
N = 100 and = 0.4(o), N = 200 and (j> = 0.2(D), 
N — 400 and <f> = 0.1(0)- Inset: Heating-dissipation ratio 
q for N — 800(o), N = 400(D), N = 200(0) and N = 100(*) 
for several values of <j). The line is a fit of the form (Nrfi) 1 ^ 2 ■ 

It has been speculated that the non-Gaussian distri- 
butions are caused by spatial correlations in the gas |13| . 
However, we can qualitatively reproduce the different dis- 
tributions we observe in simulation with a simple model 
(based on Rcf . [l4( , but including dissipation) of a two- 
dimensional inelastic gas of N particles without spatial 
degrees of freedom. In this model, every time step C pairs 
of particles are selected at random and collide using equa- 
tion ^with a random impact parameter — 2r < b < 2r, 
where r is the radius of the particles. At the same time H 
particles are selected at random to be heated by adding a 
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FIG. 5: —]n{—ln[P(v x /a x )]} versus \n(v x /a x ). The sym- 
bols shown are velocity distributions acquired by simulation 
for q = 120 (o), 0.08 (□), 0.012 (o). The lines show the ve- 
locity distributions found in the model for the same values of 
q (solid, dotted, dashed). 

random velocity. This gives us a heating-dissipation ratio 
of q = H/C. In Fig. [5] we compare velocity distributions 
from the model and simulations for different values of q. 
We find good qualitative agreement between simulation 
and model, even though it has no spatial degrees of free- 
dom. This implies that spatial correlations play a minor 
role, if any, in P{v). 

Due to the idealized nature of our system, it is not pos- 
sible to do a direct comparison between our simulation 
and experiments. Yet, some experiments seem to show a 
similar behavior, accounting for the paramet er q that we 
introduced here. In the experiments of Ref. [12j velocity 
distributions go from non-Gaussian to Gaussian when a 
rough plate or a layer of heavy particles is used instead 
of a flat plate. In this case, energy is injected directly 
into the directions parallel to the plate, effectively in- 
creasing the number of in-plane heatings over collisions. 
Most convincing is the experiment by Blair and Kudrollli 
[Tol| . Here the number of collisions is increased by adding 
more particles. As a result, their velocity distributions 



develop the same crossover we see both in our simulations 
and model. 

We studied behavior of the velocity distributions of 
granular gases as a function of </>, the area fraction, and 
77, the coefficient of restitution. Specifically, we consid- 
ered the effect of driving the gas by heating uniformly, as 
is assumed in theory and many prior simulations, and by 
heating through a boundary, as is done in most experi- 
ments. We find that there exists clear qualitative differ- 
ences between the velocity distributions for uniform and 
boundary heating, which demonstrates that the form of 
the distributions is not simply a function of material pa- 
rameters (e.g., 77). Furthermore, we show that there is no 
evidence for a universal velocity distribution with a con- 
stant exponent a = 1.5. Instead, for boundary heating, 
we find that velocity distributions cross over from one 
exponent to another for the high-velocity tail. For this 
high-velocity tail we observe a wide range of exponents 
and we find a = 1.5 only for specific values of <f> and 77. 

Instead, we demonstrate that the distribution of ve- 
locities for dissipative gases, while not universal in form, 
depends only on two parameters: the coefficient of resti- 
tution r) (a material parameter) and q = Nh/Nc, the 
average ratio of heatings and collisions in the gas (a func- 
tion of experimental conditions). We find that velocity 
distributions range from Gaussian for q ^> 1, where heat- 
ing dominates dissipation, to strongly non-Gaussian for 
5 « 1, where the dynamics of the gas is dominated by 
the dissipative collisions between particles. Furthermore, 
a simple model of a driven, inelastic gas without spatial 
degrees of freedom reproduces the entire family of veloc- 
ity distributions we find in simulation, as we vary 77 and 
q. This means that the velocity distributions are non- 
Gaussian not because of spatial correlations. Rather, it 
is the cascade of energy from a few high-energy parti- 
cles to the slow-moving bulk of the gas that is the key 
determinant of the non-Gaussian velocity distributions. 
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